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Abstract Finite-temperature diamagnetic properties of doped antiferromagnets as mod- 
eled by the two-dimensional t-J model were investigated by numerical studies 
of small model systems. Two numerical methods were used: the worldline 
quantum Monte Carlo method with a loop cluster algorithm (QMC) and the 
finite-temperature Lanczos method (FTLM), yielding consistent results. The 
diamagnetic susceptibility introduced by coupling of the magnetic field to the or- 
bital current reveals an anomalous temperature dependence, changing character 
from diamagnetic to paramagnetic at intermediate temperatures. 

The dc orbital susceptibility of the system in the external magnetic field is 
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d 2 F _ xo 



Z da 2 \Z da 
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where xo = ^oe 2 a 4 /h 2 and a = eBa 2 /h. In the previous studies [1] it was 
realized that results are quite sensitive to finite-size effects, so we also used the 
QMC method, where much larger lattices can be studied. 

The magnetic field introduced into the t-J Hamiltonian via the Peierls con- 
struction, affects only the hopping of the electrons. Within FTLM the results 
are obtained by the numerical derivation with respect to a. 

Using the standard Trotter-Suzuki decomposition and the woldline repre- 
sentation of the QMC for the fermionic models, magnetic field enters matrix 
elements concerning the hole hopping. The plaquette weights along the hole 
worldline obtain an additional phase factor. Taking the field derivatives ex- 
plicitly the expresson for the orbital susceptibility can be written as 

Xd = -Xo^ 1T l , (1.2) 
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Figure 1.1 Orbital susceptibility \d vs. T (left) for one hole obtained via QMC (dots) andFTLM 
(line) for J = OAt; different J and 7 (right). For comparison also results of high-temperature 
expansion (HTE) for J = are shown (dash-dotted line). 



where S is the projected area of the hole worldline. Xd can be thus measured 
without the presence of a magnetic field. This is just another consequence of 
the more general fluctuation-dissipation theorem. In doped systems we are 
hindered by the well known "fermionic sign problem" of the QMC, not present 
in the undoped case. Even though S 2 is strictly positive, the thermal average 
(<S 2 ) can become negative because of correlations between the Monte Carlo 
sign and the magnitude of the area S. For QMC the sizes of considered systems 
were 6x6 and 8x8. 

With FTLM a few mobile holes on a system of tilted squares with up to 
20 sites and periodic boundary conditions were considered. It is nontrivial to 
incorporate Landau phases corresponding to a homogeneous B, being at the 
same time compatible with periodic boundary conditions. This is possible only 
for quantized magnetic fields. 

In Fig. 1.1, Xd obtained via both methods is presented. For T 3> t, the 
response is diamagnetic and proportional to T~ 3 as well as essentially J- 
independent [1]. The most striking effect is that the orbital response below 
some temperature T v turns from diamagnetic to paramagnetic, consistent with 
the preliminary results obtained via the FTLM [1]. In order to locate the origin 
of this phenomenon, results for different J and anisotropics 7 are also shown. It 
appears that T p scales with 7 J, i.e. at J = the response is clearly diamagnetic 
at all T, and for 7 = 0, J > no crossing is observed with either method. 

At lower temperatures T < Td <C T p , the diamagnetic behavior is expected 
to be restored. This follows from the argument that at T — > a hole in an 
AFM should behave as a quasiparticle with a finite effective mass, exhibiting a 
cyclotron motion in B 7^ 0, leading to Xd(T — > 0) — > —00 [1]. Numerically 
it is easiest to test this conjecture for a single hole and 7 = 0. This is also true 
for J = 0. 
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Figure 1.2 \d vs - c h for several T and J = QAt. The last graph contains also a 4 th order HTE 
result (dotted). Canonical (dots) and grand-canonical (line) values for all Ch are obtained with 
FTLM on 16 sites. 



In Fig. 1.2 also results for \d for finite doping Ch > are presented. For 
nearly empty band > 0.7 the \d 1S diamagnetic and weakly dependent 
on T. In this dilute regime strong correlations are unimportant, thus Lan- 
dau diamagnetism is expected. At moderate temperatures T > J and for an 
intermediate-doping 0.2 < Ch < 0.7 the \d is dominated by a paramagnetic 
response. There is a weak diamagnetism at Ch < 0.2 and T > T p , while 
the paramagnetic regime extends to Ch = for T < T p . For low tempera- 
tures T« J quite pronounced oscillations in Xd{ch) appear and can be partly 
attributed to finite-system effects. 

The explanation can go in the direction proposed by [5], that at low doping 
c/,^0we are dealing with quasiparticles (with a diamagnetic response), being 
a bound composite of charge (holon) and spin (spinon) elementary excitations. 
The binding appears to be quite weak and thus easily destroyed by finite Torch, 
enabling the independent and apparently paramagnetic response of constituents. 
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